EPJ manuscript No. 

(will be inserted by the editor) 



OO 



o Quantum thermodynamic processes: A control theory for 



machine cycles 



^ Jan Birjukov^''*, Thomas Jahnke^, and Giinter Mahler^ 

^ Chair for Theoretical Physics and Applied Mathematics, Urals State Technical University, Mira 19, 620002 Jekaterinburg, 
Russia 

^ Institut fiir Theoretische Physik 1, Universitat Stuttgart, Pfaffenwaldring 57, D- 70550 Stuttgart, Germany 



o 

S 

> 

a 

B 

o 

(N 
> 

m 
in 
o 

O 



X 



Received: 8 November, 2007 / Revised version: 20 May, 2008 

Abstract. The minimal set of thermodynamic control parameters consists of a statistical (thermal) and 
a mechanical one. These sufhce to introduce all the pertinent thermodynamic variables; thermodynamic 
processes can then be defined as paths on this 2-dimensional control plane. Putting aside coherence we 
show that for a large class of quantum objects with discrete spectra and for the cycles considered the 
Carnot efficiency applies as a universal upper bound. In the dynamic (finite time) regime renormalized 
thermodynamic variables allow to include non-equilibrium phenomena in a systematic way. The machine 
function ceases to exist in the large speed limit; the way, in which this limit is reached, depends on the 
type of cycle considered. 

PACS. 05.30.-d Quantum statistical mechanics - 05.70. Ln Nonequilibrium and irreversible thermodynam- 



1 Introduction 

Thermodynamics [HElll] has long since been supposed to 
be an extremely efficient phenomcnological theory, though 
with regard to large systems only. However, it has recently 
been shown that thermodynamical properties emerge al- 
ready for small quantum systems provided they are em- 
bedded in some appropriate environment [4]. As a con- 
sequence, the idea that a single quantum object might 
well be described by means of thermodynamic concepts 
should no longer be considered self-contradictory. Nano- 
thermodynamics [5ll6] is becoming an emergent field that 
might be relevant to many areas ranging from nano-physics 
to molecular biology. 

Among the basic concepts of thermodynamics there 
are two standing out for their fundamental and practical 
significance: heat and work. It is not surprising that issues 
concerning quantum thermodynamic machines attract so 
much theoretical interest [TllSll^nifnifT^fnifTSlfTClfTTj. 

Various possible levels of description can mainly be 
classified by the extent to which explicit quantum effects 
are being taken into account. So, time dependence may 
either be included by means of a classical driving system 
or by coupling to another full quantum mechanical sub- 
system |17) . The dynamics of the machine may or may not 
be allowed to manifest some quantum mechanical coher- 
ence. Anyway, in this context one should observe that the 
coupling to heat baths is not an unavoidable nuisance, but 
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an essential part of the operation of any thermodynamic 
machine; decoherence |18| should thus be dominant. 

In this paper we intend to explore the universal limi- 
tations for the performance of a quantum thermodynamic 
machine as it arises from a conceivable driving scheme and 
those quantum properties of the working substance, which 
could survive strong decoherence: the discrete spectrum 
and probabilistic character of the state. For this purpose 
we employ a kind of control theory approach with two for- 
mal parameters, one related to the spectrum and another 
to the state of a quantum system (cf. piQllTO]). They are 
considered as the external control. In the ideal case this 
is all we need to capture the essentials of thermodynamic 
machines, their possible cycles and efficiencies. This sim- 
ple control model is then extended to include an internal 
time scale (phenomcnological relaxation time) to allow for 
non-equilibrium effects. This will give us a direct access 
to finite time thermodynamics [9l[T9] . 

Such an approach is entirely within the spirit of ther- 
modynamics, which itself can be formulated as a powerful 
control theory [2]. We assume that this system theoret- 
ical scheme captures the main features of any concrete 
implementation, for which coherence [8] does not play a 
major role. While there have been speculations about dra- 
matic differences between classical and quantum machin- 
ery |20j , (even claiming violations of the second law) , we 
come here to the opposite conclusion: Classical and quan- 
tum mechanical machines behave essentially identical - 
to an extent, which is almost unbelievable. In fact, this 
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scale-invariance could hardly be expected, if one took it 
for granted that the pertinent (thermodynamic) concepts 
would, indeed, only apply in the thermodynamic limit of 
the system considered. 

To be sure, we consider, in a sense, the ideal case. 
There are other important limitations in the nano-domain: 
as all length scales shrink, also the possibilities of ther- 
mal isolation become severely constrained (cf. [21 ). This 
means, e.g., that the interaction with baths of different 
temperatures may no longer be assumed to be switched 
on and off at will. Leakage becomes unavoidable 11,12| 
[TT] . This aspect will presently be excluded, as are any 
implementation issues. 



2 Model of Control 



at all times, i.e. the so-called quasi-static limit [7. 

In the very general case, all Ef'^'^{'^) had to be regarded 
as N independent functions that would make the model 
hardly tractable. In what follows we consider, instead, the 
special class of the spectral control 

Eril)^9il)-e., (3) 

where g is some monotonous function independent of i, 
and {£i}fLi are to be regarded as a set of characteristic 
constants. Two pertinent examples resulting in 17(7) = 
7 and (7(7) = 7~^ are given in Sec 13.5.11 and Sec 13.5.21 
This control allows to treat the distribution pi (a) formally 
independent of 7 even if the spectrum does explicitly enter 
the occupation numbers as, e.g., in the canonical case, 

cf. gsi). 



We consider a quantum system with discrete spectrum 
embedded in some environment. One can imagine three 
qualitatively different functions the environment could per- 
form with respect to the system: a mechanical control, a 
decohering bath and a thermal bath considered clS Si stS'" 
tistical control. 

Provided the weak coupling conditions are met, we as- 
sume that the mechanical control comes into action through 
a parameter-dependent spectrum of the effective Hamilto- 
nian, {£'f**(7)}^i, where N is the number of levels. One 
may regard this control as purely mechanical as long as 
the level occupations can survive the corresponding spec- 
trum disturbance. It proves to be the case in the presence 
of decoherence rapid enough (time scale r^ec) compared 
with 7 dj. 

This kind of decoherence can well originate from the 
appropriate part of the interaction between the system 
and its environment. In our model we assume such a de- 
cohering bath as a physical prerequisite of an adiabatic 
process. It has another important consequence, effectively 
making the coarse-grained density matrix pt and the mo- 
mentary hamiltonian _ff®**(7i) commute, excluding as well 
autonomous Schrodinger-dynamics. So we have to assume 
that our quantum system always stays in a mixed state 
determined by the diagonal elements of p in the respective 
energy representation, {pi = Piilfli. 

In order to introduce the statistical control, we postu- 
late the existence of a one-parameter family of distribu- 
tions, {'Pi{a)}iLi: inherent to particular conditions of the 
contact between the system and its environment. These 
distributions are assumed to be stable in the sense that 
whenever the actual distribution pi differs from Pi{a) at 
some given a, this deviation will decay on a time scale 
tr » Tdec according to 



Pi 



(1) 



(relaxation time approximation). We shall refer to such 
Pi{a) as attractor. 

Note that for small compared with the characteris- 
tic time of enforced a-parameter alterations, one gets 

Pt=Pt{a) (2) 



3 Quasi-static limit 

The quasi-static limit plays a fundamental role in stan- 
dard thermodynamics. Here it describes an important ref- 
erence scenario: the limit of perfect control over a quantum 
system. Both the spectrum ([3|) and the energy distribu- 
tion ([2) are completely specified by 7 and a taken to be 
external control parameters. Hence, as long as the thermo- 
dynamic quantities are properly defined, one can consider 
an arbitrary point in the (a, 7)-plane as a thermodynamic 
state implemented on the system under consideration. 



3.1 Thermodynamic quantities 

As follows from the foregoing, for the system's state, cor- 
responding in the quasi-static limit to a certain point in 
the (a,7)-plane, one can define the entropy S and the 
internal energy U [3] : 



S{a) := - '^Inpiia) -piia), 
U{a,-f) := gij)^e^ ■ p^{a); 



(4) 
(5) 



From ((4]) we observe, taking into account the normal- 
ization condition for {pi{ct)}fLi, that 



dS 
da 



dpt 
da 



(6) 



In order to introduce the notion of temperature, we 
employ its formal definition ^i3j as the conjugate variable 
to the entropy S: 



T := 



dU 
dS 



dU\ fdS 



da J ^ \da 



Together with IH and ^ one gets 

T{a,-f) = g{j)0{a)-\ 
where 0{a) stands for 

X;inpi(dpj/da) 



0{a) := 



X;Q(dpi/da) 



(7) 

(8) 
(9) 
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a function depending solely on a. 

Consider now the total differential of the internal en- 
ergy Is]) as a function of these control parameters [7]: 

d^(-'7)=(f)^d.+ (^)^d,. (10) 

Based on ([7]) we introduce the infinitesimal increment of 
heat, dQ = TdS: 

As to the second term in pUj) . its meaning is the internal 
energy increment under constant a, in fact under constant 
entropy, cf. ([4]). It allows us to consider this term as work: 

3W{a,j) -.^ (j^^ d7=^^e,;p,d7. (12) 

Here, all the thermodynamic quantities have been de- 
fined based on a specific process in the control plane. Typ- 
ical measurements (of temperature T, say) would rather 
exploit the dependence of a specific observable A, like the 
magnetization of a paramagnetic salt, on T. 

3.2 Specific processes 

So far as every point in the control plane (01,7) repre- 
sents some thermodynamic state, a continuous sequence 
of states, a process, can be defined by a path in this plane, 
eventually by some constraint on the accessible values of 
control parameter pairs, /(a, 7) — const. For a concrete 
implementation the class of actually feasible processes will 
severely be restricted, just like for macroscopic machines. 

Within our control theory framework the simplest ex- 
amples for ideal conventional processes are isentropes and 
isochores, identified with a = const and 7 = const, re- 
spectively. It is noteworthy that along an isentrope the 
temperature ^ is directly proportional to g{j). With 
(7(7) = 1/7, for example, T(a,7) will decrease with in- 
creasing 7: 

Tex- (13) 

7 

This is what happens to photons in a cavity of size 7 = L, 
a phenomenon known also from the photon temperature 
reduction in our expanding universe (cosmic microwave 
background [22] ) . Another example is adiabatic demagne- 
tization (magnetic cooling), with j — B and 3(7) = 7 (see 
Sec 13. 5. 1| resulting in T c>c 7. 

Based on ([8]) , the definition of an isothermal process is 
straight forward: T{a,"f) = const (~ T). This constraint 
can be cast into the convenient form: 

5(7)=T.e(a), (14) 

In which T plays the role of a scale factor: The shape of 
every isotherm is the same, depending eventually on the 
chosen form for 5(7) and {'Pi{c()}^i, the attractor. 



If, as a particular choice, one takes the canonical at- 
tractor: 

P.(a) = ^Ca„e-"^% ^Ca„ = ^e-"^'; (15) 
one finds from ([9]): 

Ocania) = a. (16) 

With 3(7) = 7 the isotherms \1A\ will now be just straight 
lines in the control plane, 7 = aT at constant T. With 
5(7) — 7^^ would get 7 = 1/ \/aT, a sort of hyperbolae 
(see Fig. Hand Fig.[2|. 

3.3 Heat and work 

Bearing in mind the forthcoming analysis of various ther- 
modynamic cycles, which underlie heat engine or heat 
pump operation, we consider the heat and work along spe- 
cific processes. Integrating pT]) and ([T^ . one gets for the 
isotherms 3(7) = T • 0{a): 

Qr = T^e, j 0{a)dp,{a)- (17a) 

/A 
P^{a)^da■ (17b) 

for the isentropes a — const; 

Qs = 0; (18a) 

Ws = Ag{^)Y,e,p,(a)- (18b) 
and for the isochores 7 = const: 

Q7 -.9(7)5I«»^P'("); (19a) 

= 0; (19b) 

where A denotes the corresponding increment along the 
process line. 

The sign of the heat flows Qr and is calculated 
in Appendix [B] for the case of the canonical Pi{a) p^]) . 
Such general sign statements cannot be proven for the 
work inputs Wt-s, as their sign depends, together with 
the rest, on the particular choice of {eilfl^. 

3.4 Carnot cycle 

Provided the RHS of ([9]) behaves well, it is always possible 
to compose a closed path in the (a, 7)-plane from two isen- 
tropes and two isotherms - thus leading to a conventional 
Carnot cycle. 

All the pertinent thermodynamic quantities, such as 
entropy (H]), internal energy ([5]), temperature ([5]), heat pT]) 
and work (|12p are defined in the quasi-static limit in such 
a way that Eq. (flO]) turns into the Gibbsian fundamental 
form: 

AU ^TAS + (1W, (20) 
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Fig. 1. Carnot's cycle in the (a,7)-plane associated with the 
canonical attractor (jlSp and 3(7) = 7. It is composed by 
segments of isentropes a — 0:1,2 and isotherms 7 = QTc.h; 
(Th > Tc). Arrows correspond to heat-engine performance. 
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Fig. 2. As Fig. [T] but for (7(7) — 7"^. The Carnot cycle is 
composed by segments of isentropes a — ai.2 and isotherms 
7 ~ 1 /y^aTcy, {Th > Tc). Arrows correspond to heat-engine 
performance. 



irrespective of 3(7) or the chosen kind of attractor ([2]). 

Thermodynamic efficiency of a heat engine cycle is de- 
fined as the ratio Wont/Q%, where Wout = —Wo > is 
the work output per cycle and Qti, > is the heat input 
during the isothermal stage at the higher of two assigned 
temperatures: Th > Tc. As usual, Qr^ < is supposed 
to be discarded, i.e. cannot be re-used. The validity of 
the Gibbsian fundamental form immediately leads to the 
Carnot efficiency [2] 



1 



Tc 
Th' 



(21) 



for any heat-engine cycle of Carnot's kind in the {a, 7)- 
plane independent of such model details as {pi{a)} or {ei} 
or g(7). 

Thus, the Carnot efficiency as a limiting fundamen- 
tal value does not even require thermal states, only the 
control scheme with two parameters, as stated. The tem- 
peratures Tc,Th would then only have a formal mean- 
ing, though. In some sense, this is a generalization of 
the well-known universality established in conventional 
macroscopic thermodynamics. 

In order to present an illustrative example, we explore 
the case of the canonical attractor (jlSp . The corresponding 
Carnot cycles are sketched in Fig. [Hand Fig.O 

According to ([SO]) , for the heat-engine performance one 
must drive the cycle in Fig.[l]anticlockwise, while the cycle 
in Fig.[2]clockwise: These directions correspond to heat in- 
put for the high- and heat output for the low-temperature 
isothermal stage, respectively. Treating both cycles step- 



wise based on (fT7|) and (fT8|) , one explicitly gets the Carnot 
efficiency (|2T|) . as expected. 

Two examples for other possible choices of the attrac- 
tor type are given in Appendix |^ These examples illus- 
trate the diversity of forms the Carnot cycles can have 
in the control space (a, 7) depending on the underlying 
attractor. In the quasi-static limit, however, all of them 
share the standard rectangular form in the (S*, r)-space 
(Fig. 2]). Also the efficiency is always given by the Carnot 
value; in this sense, one cannot win anything by trying to 
implement exotic distribution functions. 



3.5 Otto cycle 

It is well-known that the universality established for the 
Carnot cycle does not carry over to other cycle types. This 
means that the resulting efficiencies rj < rf would then 
depend on the actual control functions 5(7) and pi{a). 

The Otto cycle is in some sense the most fundamental 
quantum thermodynamic cycle: On the isentropes, there 
is only a change of the spectrum while on the isochors 
only the occupation numbers are changing. Most of the 
theoretically discussed quantum thermodynamic machines 
perform, in effect, the Otto cycle [8lfT5ll23j . 

Fig. [3]shows the Otto cycle in the (a,7)-plane. As the 
Carnot cycle discussed above, this cycle must be driven 
anticlockwise or clockwise, depending on the choice of 
(7(7), in order to get the heat-engine performance. 

If the function 3(7) is increasing, ([18]) and ([T9|l lead to 
the efficiency: 

5(71) (22) 



= 1 



5(72) ' 



For a 2-level system, one gets, identifying 

^{1) = 5(7)(e2 - ei) 
the well-known result |24| 



^(71) 
^(72)' 



(23) 



(24) 



which is often mistaken as representing the Carnot effi- 
ciency [25] . 

With decreasing 5(7) one must drive the cycle clock- 
wise to get a heat-engine. The efficiency is then given by 



gill) ' 



(25) 



In any case, the efficiency only depends on the relative 
compression, just like for the classical Otto cycle [1]. 

To compare these results with the efficiency of the 
Carnot cycle, we consider the canonical case ([m[TB]) to 
get 

o -, Tca2 



77- 



1 



Thai 
Th^i 



(26a) 
(26b) 
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Ql 



Fig. 3. Otto cycle in the {a, 7)-plane. The isentropes are given 
by Of = const and the isochors by 7 = const. The dashed lines 
are canonical isotherms with the highest and lowest tempera- 
ture of the cycle for (7(7) = 7, while the dotted ones hold for 
9(7) = 7"^- 



where Tc (T'^) is the lowest temperature and Th (T^) the 
highest temperature along the cycle. As shown in Fig. [3l 
Tc is reached at point 2 and Th at point 4 for increasing 
3(7), whereas for decreasing 3(7) Tj, is reached at point 3 
and Th at point 1 (cf. f^). As one can see, the efhciency of 
the Otto cycle remains below the Carnot value. Formally, 
the Carnot efficiency is reached for a2 = ai, when the 
total work output of the machine vanishes (cf. [15]). 

In the following we shall give two examples of spectral 
control: 



4 Driven non-equilibrium 

Now we are going to re-define the thermodynamic quan- 
tities introduced in Sec. [31 allowing for deviations of the 
momentary distribution {pi} from the attractor {pi{a)}. 
In this way we extend the scope of our consideration be- 
yond the quasi-static limit ^ including processes with 
the cycle times comparable with the relaxation time tr 
introduced in P^. 

4.1 Beyond the quasi-static limit 

According to ((T|), the probabilities pi no longer depend 
just on the momentary control a(t) as in the quasi-static 
limit; they also become dependent on their past history. 
Nevertheless, they are quite well computable, as soon as 
a{t) and the initial values are defined, and one can still de- 
fine the entropy and the internal energy just by replacing 
Pi{a) in ([U and ([5]) with pi. 

We define the non-equilibrium work dW* as the in- 
crement of the internal energy due to variations of the 
mechanical control 7: 

while the non-equilibrium heat SQ* based on the first law: 
dQ* =dU* -aw* ^giJ)Y,^^<^P^■ (32) 



3.5.1 Example: Spin in a magnetic field 

For a spin in the magnetic field B, Ei (x B holds. Identi- 
fying the control parameter 7 with B, we have 



5(7) = 7 

and therefore the efficiency (|22p simply reads 



r7°- 

'/spi 



1-^, 
72 



(27) 



(28) 



3.5.2 Example: Particle in a box 

Consider a particle in a box of length L [20] . The energy 
levels depend on the length as Ei oc 1/L^. This length can 
thus be used to control the spectrum. Identifying 7 with 
L, we have: 



= ^ 



(29) 



Because 3(7) is decreasing here, we have to take ([25]) as 
the efficiency: 



Vho 



(30) 



(By an asterisk we indicate the renormalization due to 
non-equilibrium. ) 

The renormalized temperature T* will be taken to re- 
main a control variable asking for the response of U* to 
a change of the statistical control a. Note that the non- 
equilibrium distribution {pi} can formally be represented 
as 

Pi{a, Api) = pi{a) + Api, (33) 

with deviations Api due to the dynamical response of the 
system. Being well computable by means of the relaxation 
equations ([1]), they are, however, not under our direct con- 
trol. That is why we seek the response to a change of a at 
some given {Api}: 



da 



1 i^Pi 



dS* 
da 



-1 



(34) 
(35) 

The quantity T*, determined in this way, we shall call 
renormalized process temperature. 

In the following we assume the time-dependence of the 
control parameter a{t) to be: 



Such a definition results in (cf. 

J2<^i{dpi/da) 



X;inpi(dpi/da) ' 



a = aQ + vt, 



(36) 
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where v may be positive or negative. Of course, there are 
other possible choices, indicating that more details of the 
control may become important. 

With the attractor {pi{a)} driven according to ([36|l . 
the relaxation equations ([T]) take the form: 

(dp^/da) = -{vTj^y^{p, -pi{a)) (37) 

and yield the solutions: 

=p(°)e-^ + (uT^)-M p,{a')e~'^Aa'. (38) 



This means that, in a sense, the pi could still be considered 
as an explicit function of a as soon as their initial values 
p1 and ao are specified by the given process history. 

4.2 Renormalized process temperature 

In order to understand what happens with the process 
temperature along non-isentropic paths in the control plane, 
we consider the asymptote under very slow driving, 
{vTfiJAa) ^ 1, when the following expansion holds (see 
App.P: 



p^ =Pt[a) -'"'^R-^ + 



Aa 



(39) 



where Aa denotes the increment of a in the course of the 
process. Inserting this expansion into psp and keeping 
terms up to first order in {vTj^/Aa), one gets: 



T* 



T{a,l) 1 + 



VT 



E(l/K)(dp»/d«J 
« E In p, (dp, /da) 



Aa 



(40) 

observing as well the definition ([5]) for the temperature 
T(a,7) in the quasi-static limit. 

Of course, any chosen process line f{a,"f) — const 
implies a certain relationship between 7 and a. In the 
canonical case (IT5|) and when the process is taken to run 
along the line (7(7) = T • a, which would be an isotherm in 
the quasi-static limit, cf. ([Till . Eq. ((40|) reduces to: 



(41) 



T*{a) = T 1 + — ^ 



a 



Aa 



For a Carnot machine cycle (Fig. [T] and Fig. [2]) this 
means, compared with the quasi-static limit, that the dy- 
namically renormalized process temperature is increased 
for increasing a (v > 0, contact with a heat sink) and de- 
creased for decreasing a (u < 0, heat source). Temperature 
shifts like those following from (|¥T|l have ad hoc been in- 
troduced in Curzon and Ahlborn's finite time analysis of 
the heat engine efficiency |19j. Already at vtr « 1, how- 
ever, the ST loop proper is dominated by counteracting 
non-equilibrium excursions (Fig.lH cf. [9^), which indicates 
that there will be no machine function left at high driving 
speed (cf. Sec lS.ip . 

One may wonder, if such a behaviour of the tempera- 
ture as shown in Fig. 3] really can be observed in physical 
systems. Indeed this seems to be the case: Very recently 
such a temperature pattern was found in numerical exper- 
iments for the finite-time Carnot cycle of a classical ideal 
gas by means of molecular dynamics simulations [27]. 



. 50 
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1 3* 



. 72 



1.25 



Fig. 4. Carnot cycle in the (5**, r*)-plane. Both cycles from 
Fig. [1] and Fig. [2] look similar here. The rectangle denotes the 
perfect Carnot cycle in the quasi-static limit vtr = 0. The solid 
line corresponds to a stationary cyclic regime at vt^ = 0.09; 
the dashed one - at tjt^ = 0.9. This is a result of numerical 
simulations in accordance with S* :— — J^lnpi • Pi and psp . 
with Pi under relaxation 
chosen here to be {—3/2; 



08}. The set of constants {ei}fLi is 
-1/2; 1/2; 3/2}. 



4.3 Renormalized heat and work 

Consider the renormalized heat and work along those spe- 
cific processes underlying Carnot and Otto cycles. One 
has to integrate now ([32)1 and pT|) along the correspond- 
ing lines in the [a, 7)-plane. In the case of the canonical 
isotherm 5(7) — aT this leads to the integrals J^^ a'dpi 

and J^^pida', with the probabilities pi under relaxation 
as given by ([55)1 . The latter is readily achievable directly 
from the relaxation equation (|37p : 

/ da' = I p,{a')da' - vtj^ U ~ pf) , (42) 

immediately followed by: 

i pa 

a' dpi = / a'dpi{a') + a {pi - Pt{a)) (43) 

Jao 

- "0 (pf ^ - Pj(ao)) + vrj^ (p, - pf 



Renormalized isothermal heat and work are thus in the 
canonical case: 



Qt = Qr + q + q' + q\ 

q' := -aoT^e, (^pf^ - p,{ao)j , 
q" := aT^e, (p., - Pt{a)) ; 



Wr ^ Wt + '> 
w := —VTfiT 



Pt - Pi 



(o)\ _ 



(44a) 



(44b) 



Again, the respective renormalized quantities are denoted 
by an asterisk, while the heat and work in the quasi-static 
limit (right hand side) are defined according to (fTT]) with 
0{a) = a. 

Integration along the isentropic line a = const yields 
zero for the heat and Ag ^ ^ipf^^ for the work, due to 
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the conservation (ideal in our model) of the respective 
distribution — the results arc similar to those obtained in 
the quasi- static limit (fTS)) : 

Qs = 0; (45a) 
Ws = Ws+w', 

w' ■.= Agi^)Y,e^{p["^-Mao)); (45b) 
Finally, for the isochores 7 = const one gets: 

Q* = Q-, + q + q, 

q 9{l)^ei{Pz-Pi((^))- (46a) 
W; = 0; (46b) 

where Q*, as well as Wg before, is formally expressed 
through its quasi-static value. 



5 Finite-time thermodynamic cycles 

In order to realize a cyclic process in the (a, 7)-control 
plane, we have to assign a{t) for both running directions 
between the turning points ai and a2- We allow for dif- 
ferent driving speeds, namely: 

ai-,2{t) ^ ai + nt; a2^i{t) ^ a2 ~ Xnt, (47) 

with K, A > 0. 

Taking ([38)l with v = k and v = — Ak, respectively, one 
comes up with the following distributions at the turning 
points in the stationary cyclic regime: 



Pi — (e^"^« — e 1 • (kt^)^ 

Pi(a') ( e _ e ^"^H ) da'; 



(48a) 



.(2) 



(48b) 



da'; 



5.1 Finite-time Carnot cycle 

Consider now the Carnot cycle (Fig. [1] and Fig. [2]) driven 
with the speed t; = k on the low- and v = — Ak on 
the high-temperature isothermal stage according to (|T7|) . 
The established cyclic regime implies here the distribu- 
tion (|i5a|) at the points "1" and "4" and (|48bp at the 
points "2" and "3" of the cycle — there is no relaxation 
on the isentropes, as discussed before. 




Fig. 5. Corrections (149 [1 and (|50p to the quasi-static heat and 
work along the stages of Carnot- like machine cycle (Fig. |4]). 
Straight lines denote work, wavy ones - heat. Arrows pointing 
into the contour indicate energy flows into the system. (The 
primary, quasi-static contributions (|17p and (|18p providing the 
machine action are not shown.) 



With (|44|) we get the following correction terms to the 
quasi-static heat and work along the isotherms 1 — )■ 2 and 
3 — > 4; respectively: 



912 


= '«Ti?Tc X ei {pf^ 




9'l2 


= -aiTc X e, (pf ^ 


-Pi(ai)) > 0; 


9;'2 


= a2Tc X ei [pi^ - 


p,(a2)) > 0; 


934 


= A(Th/Tc) 912 < 0; 




934 


= -(VTe)gi'2 <0; 




9^'4 


= -(Th/Te)gi2 <0; 




W12 


= -912 > 0; 




W34 


= -934 > 0. 





(49) 



Eq. (j45b[l yields the correction to the quasi-static work 
along the isentropes 2^3 and 4^1, respectively: 



^3 = a2ZiTXe, (pf^ -p.(a2)) > 0; 

Z\TXe,(p,«-p,(ai)) >0; (50) 



-ail 



Here we have observed Zi.9(7) = aZ\T, which follows from ([14] 
and (fTe]) . 

The inequalities (sign) in (j49p and (|50p are proven in 
Appendix |B] The sign of every correction term remains 
the same for any k, and A, allowing us to sketch their 
general scheme in Fig. [SJ Of course, all the terms listed 
here relate to the corresponding cycle step as a whole, 
not allowing for physical separation. Nevertheless, such a 
representation proves to be useful for the purpose of the 
heat and work transfer analysis in various driving speed 
regimes as well as in comparison with the cycles of other 
kinds. 

In the limit of very slow driving, (KT^/Z\a) ^ 1, one 
has (see App. [C|): 



pf^ -P^P = Pi{a2) - Pi{ai) +0 [kt^] , 



'R\ 1 



(51) 
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which means that all the terms in (|49|) and ([50| are of first 



order in kt^. In Sec. l5.2l these results will be used in order 
to reveal the fast degradation of the machine efficiency 
with increasing driving speed k. 

In the limit of very fast driving, {Aa/nTj^) <C 1, the 
terms above can be shown to be (App. [C|) : 



(2) (1) 
Pi -Pi 



(i;2) ~ ( \ 
Pi -Pi{ai;2) 



Aa 



1 - 



(52) 



Pi -Pi(ai;2) + 



Aa 



where 



2 ro'2 

li := — — / (pI-pi(a))Q!da, 



1 
Aa 



Pi(a)da; 



(53) 



- the latter being the average of Pi{a) over the interval 
[ai; 02]. 

The result (|52p supports the intuitive expectation for 
the collapse of the cycle on (5*, T*)-space (Fig.|3]) at high 
driving speed k. The only energy flows surviving kt^ 00 
are the total work along the isentropic stages, W2*3-4i; ^^'^ 
the quasi-static part of the work along the isotherms: 



W^2*3;4i = ± ^Tq;2;i ^ e^p,; 

M^12;34 = ± Tc;h^a ^ ^iPl- 

The heat flows Q* in the limit {Aa/nTj^) ^ 1 are: 



(54) 



Q12 — 



Aa 

KTr 



-Tc— y 4 

R 

Aa 



•r 



(55) 



negative on the both isothermal stages (App. |B]). Sure, 
there is no machine action left at kt^ 3> 1. 

A comment must be made concerning the scheme in 
Fig. m The triples (q'12; W41; q34) and (qia; w^g; q'34) look 
like effective heat-pumps, counteracting the primary ma- 
chine action along the cycle. This impression is supported 
by the balance relations, following from and (I5D1) : 



q'i2 + W41 + q34 = 0; 

qi2 + W23 + q34 = 0' 



(56) 



and may insinuate the idea of a crossover to a heat pump 
at some values of k and A. As a matter of fact, this never 
happens for a cycle of the kind under consideration; it is 
proven in App. |D] that QJ2 < always holds — there is 
no heat absorbtion from the cold bath as a result of the 
low-temperature isothermal step of the cycle. 

We conclude this paragraph with graphs (Fig. [6] and 
Fig. E]), which illustrate the degradation of the Carnot 
cycle's machine action with the growing driving velocity. 



^34 
1 w 



Q'l2 



(a) 



(b) 




Fig. 6. The development of the Carnot cycle can be divided 
into five steps: For small ktr, the cycle works as a heat engine 
(a). At a certain velocity, the total work vanishes, and there 
is only heat flowing from the hot to the cold bath (b). After 
that point, work has changed its sign (c). Then there exists a 
velocity, where (534 vanishes. The work therefore is completely 
transformed into heat, flowing into the cold bath (d). Finally, 
heat is flowing in the cold as well as in the hot bath (e), which 
corresponds to the high driving speed limit, discussed before. 
Note that the size of the arrows is rescaled at each picture: All 
components of heat and work tend to decrease for increasing 
velocity as shown in Fig. [T] 



0.6 - 
0.5 -' 
0.4 - 
0.3 - 
0.2 - 
0.1 - 


-0.1 
-0.2 
-0.3 




Fig. 7. Behaviour of Wo (solid lines), Q34 (dashed lines) 
and Q12 (dotted lines) with increasing velocity, typical for 
a Carnot cycle. The chosen parameters are {ei}fLi — 
{-3/2; -1/2; 1/2; 3/2} and ai = 0.2, Q2 = 0.8, Tc = 1, Ti = 2. 
As one can see, both Wo and Q34 change their sign at certain 
finite velocities yet, while Q12 always stays negative. 



5.2 The efficiency at maximum power output 

In the following we restrict ourselves to corrections to the 
quasi-static heat and work along the cycle stages linear in 
the driving speed. In the appropriate limit (j5ip eqs. (j49F 
[5(1)) result in 



{ W12; W34 } = { -qi2; -q34 } = { Tc; ATh } X ktrq, 

(57a) 



where a := {Pijai) ~ Pi{a2)) > 0; 



{ q'12; W41; q34 } = { Tc; /iT; -T^ } x XnTRb' , 

where 6' := ai ^ ei(— dpi/da)^^ > 0; (57b) 

{ qi2; W23; q34 } ^ { Tc; at; -t^ } x uTRb", 

where 6" := a2 ^ ei(— dpi/da)^^ > 0. (57c) 
(The inequalities are proven in Appendix iBl) 
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We assume, in addition, that the relative part of the 
time per cycle spent on isothermal steps remains always 
the same, say 6. Hence, the period At for one full cycle 
turns out to be [19] : 



At 



Z\a A + 1 
6k a 



(58) 



The ratio of the total work output per full cycle, W*^^ = 
—W* > 0, to the heat absorbed on the high-temperature 
isothermal step, Q'^^, yields the efficiency 77*. As follows 
from the scheme on Fig. [5] and ([57)) . one gets: 



V 



Wo^t - KT^ ((Tc + ATh)a + XATb' + AU") 
Qt, - KT^Th (A(a + h') + h") 



(59) 



One easily convinces oneself that this result is bounded 
from above by the quasi-static Carnot limit 



(60) 



reached at k = 0. The power output V* = W*^^/At is 
given by 



^kA 



(61) 



Z\a(A + 1) 

X (vFout - '«Tfl,((Tc + ATh)a + XAlb' + Alh"] 

which confirms that V* — for k = 0, i.e. for the maxi- 
mum Carnot efficiency rf . 

Our aim now is to maximize V* with respect to pa- 
rameters K, and A and to find the corresponding thermo- 
dynamic efficiency (j59p . 

The condition for an extremum, d^V* = d\P* = 0, 
yields 

Wout = 2kt^((Tc + ATh)a + XAlb' + Alb"), 
Wout = '«T^,((Tc + ATh)a + XAW + Alb") (62) 
+ KTjjX{X + l)(Tha + Alb'). 



It follows from ((62)) that the maximum of P* occurs at 
2 _ Tc l + AW'/T^a 



X' 



Th 1 + Z\T6'/Tha 



(63) 



Eliminating now Q%, Wout and A in (|59[) by means 
of (|60p and those maximum conditions, one obtains: 



VmaxV — V 1 + 




As one can see, while the formal upper bound for rj^^^j, 
is still Carnot's rj'^ , the lower bound is the celebrated 
Curzon-Ahlborn efficiency p]9l[26] : 



,7^-^ = 1 



(65) 




Fig. 8. Power output V* = W^M/^t (arbitrary units) vs. 
Efficiency ri* = Wom/Qt^- The solid line and the lower utr- 
scale correspond to the linear driving speed approximation (|59p 
and (|6ip . at A = 0.756 (|63[) . The dashed line and the upper 
KT_R-scale correspond to exact calculations based on (|48p , when 
the maximal power output appears to be achieved at A = 1.040. 



attainable in the case, when Alb' /T^^a and Alb" /Tea can 
be neglected. 

Formally, it is the triples (qi2; W41; q34) and (q'/a; w'23; q34) 
in the finite-time corrections to the heat and work (I57p . 
that makes rj^g^^-p > rf^^. In fact, they originate from 
the discrepancy between the steered attractor {pi{a)} and 
the values of pi lagging behind at the cycle turning 
points. It is clear, that such a residual non- equilibrium 
effect could not emerge in the purely phenomenological 
thermodynamic setup by Curzon and Ahlborn [19] . Nev- 
ertheless, it has been observed for the classical ideal gas 
in recent numerical experiments (27 \ . 

To be sure, the higher efficiency rj^axV does not mean 
any practical gain in the machine action here, because the 
corrections mentioned above do reduce both the heat Qt^ 
absorbed from the high-temperature bath and the work 
output W*ut per cycle, cf. (|55|) . 

Applying Taylor's formula for Pi(a) in (|57p provides 
a rough estimate for the terms Alb' /ih^ and Alb" /icO,- 
Observing and (|14p as well, one gets: 



Alb' _ {Ag)s, 



(A9) 



ATb" 



52 



(66) 



where {Ag)-r^ t, and {Ag)si 2 stand for increments of func- 
tion (7(7) along the corresponding isotherms and isentropes. 

Thus, one can expect Curzon- Ahlborn's result for 
Carnot cycles with small enough Al and large enough Aa 
(eventually large enough AS) . With the cycle parameters 
we have taken for our illustrations (Fig. [T] and Fig. ^ it 
is obviously not the case, however, and the optimum 77, 
indeed, lies between the Carnot and the Curzon-Ahlborn 
bounds, both for the exact and the linearized calculation 
(Fig.©. 



5.3 Finite-time Otto cycle 

We finally turn to the Otto machine cycle shown in Fig.O 
driven now with finite speed. Again we allow for two differ- 
ent speeds on the isochors, in accordance with ([47]), and 
no bath coupling and no relaxation are assumed along 
the isentropes. The distributions pi at the cycle's turning 
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0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Fig. 9. Behaviour of Wo (solid lines), (334 (dashed lines) and 
Q12 (dotted lines) with increasing velocity for an Otto cycle 
of a four level system. The chosen parameters are {eijfLi — 
{-3/2; -1/2; 1/2; 3/2} and ai = 0.2, = 0.8, 71 = 1, 72 = 3. 
As one can see, the ratio —Wo/Qzi, which defines the effi- 
ciency, always stays the same. 



poin ts "1 " ("4") and "2" ("3") are again given by (|i5a|l 
and (j48bp . respectively. 

The heat and work along the cycle steps then follow 
from (gSl) and (gS]): 



Q*i2 = 9{li)^^i [pl 

1^2*3 = [5(72) -.g(7i)]E^^pf^ 



(67) 



the inequalities are proven in App. [B) 

The efficiency can now easily be calculated: With in- 
creasing (7(7) we get 



23 



41) 



Q 



= 1 - 



34 



g(7i) 
.9(72) 



(68) 



With decreasing 5(7), one has to run the cycle clockwise 
to achieve heat-engine performance, which results in the 
efficiency 

.9(72) 



V 



1*0 



= 1 - 



5(71)' 



(69) 



i.e. exactly the same as in the quasi-static limit (|25p . So, 
the efficiency of the Otto cycle is independent of the driv- 
ing speed. In particular, the efficiency does not depend on 
the power of the engine, what is totally different from the 
behaviour of the Carnot cycle discussed before. The con- 
stant efficiency does not mean, of course, that the heat 
and work are constant too. All energy flows go to zero 
when the driving speed grows, as shown in Fig. [SI 

Fig. [TUl shows a ST-diagram for an Otto cycle. In this 
diagram, one can also see the decrease of the work output 
per cycle for increasing speed, since the encircled area de- 
creases. But in contrast to the Carnot cycle (Fig. d]), the 
finite speed causes spikes rather than loops. 



2 - 

1.8 - 

r6 - 

1.4 - 

1.2 - 

1 - 

0.8 - 

0.6 - 

0.4 ^ 




0.45 0.5 0.55 0.6 0.65 0.7 0.75 O.S 



0.85 0.9 0.95 



Fig. 10. ST-diagram for an Otto cycle of a four level system 
with = {-3/2; -1/2; 1/2; 3/2} and qi = 1, aa = 2, 

71 = 1, 72 = 2. The dotted lines hold for the quasi-static limit, 
the solid lines for kth = 0.5. The points 1 to 4 correspond to 
the respective points in Fig. |3l 



6 Conclusions 

We have considered a single quantum object with a dis- 
crete spectrum, an open system. The impact of its environ- 
ment has been supposed to be reducible to a parametrized 
distortion of the spectrum and to a parametrized change 
of the occupation numbers. The internal energy U and the 
von Neumann entropy 5* can then be expressed as unique 
functions of these controls. The same holds for any perti- 
nent thermodynamic quantity including those connected 
not just with a state, but with a process, i.e. heat and 
work. While there can be no operators underlying these 
quantities, they can, nevertheless, be defined for any ap- 
propriately embedded quantum system, even down to a 
single spin. 

Thermodynamic machines arise, if a cycle is enforced 
on such a two-dimensional control space. The Carnot limit 
for their efficiencies has been proven and interpreted to re- 
sult from the interplay between mechanical and statistical 
control rather than from the spectrum details or thermal 
equilibrium as such. Consequently, there is no way to vio- 
late the second law of thermodynamics, not even by using 
exotic attractor states. 

While coherence has been excluded here, dynamical 
effects have been included after defining an internal relax- 
ation time scale. The non-equilibrium features resulting 
from the finite control speed have been incorporated in 
terms of renormalized thermodynamic quantities. As has 
been shown, under such a non-quasi-static condition even 
the concept of a process temperature can still be applica- 
ble. Compared with the quasi-static one, this renormalized 
temperature turns out to be lower or higher, depending on 
the heat flow direction along the process. 

The finite time heat and work exchange as well as the 
thermodynamic efficiency have been examined in detail, 
analytically and numerically, for two standard machine 
cycles: Carnot and Otto. Though there are some qualita- 
tive distinctions, both display rapid degradation of their 
machine function, when the driving speed is increasing. 
Contrary to the Carnot cycle, the Otto cycle's efficiency 
remains speed-independent while energy flows decrease 
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Th ^ 




















a 



ai 



a.2 



1 -N"^ 



Fig. 11. Carnot cycle on the (Q,7)-plane associated with the 
nearly-uniform distribution (|70[) and (7(7) = 7. It is composed 
by segments of isentropes a — ai,2 and isotherms 7 — Tc,h • 
6>w(q) dZl); (Th >Tc). 



up to zero. For the Carnot cycle the celebrated Curzon- 
Ahlborn result proves to be the lower bound for the effi- 
ciency at maximum power output. 

We have as yet considered a special class of time de- 
pendence for the controls a{t) and 7(t), namely, a lin- 
ear one. Our approach is meant to be some kind of a 
minimal model capturing the essentials of any quantum 
thermodynamic machine deprived of coherence. Of course, 
this is not a substitute for investigations of particularized 
models; nonetheless, the approach presented allows one 
to treat any type of cycle under arbitrary dynamic regime 
on equal footing, putting aside, though, by which physical 
means this control might be implemented. 



A Carnot cycle for non-canonical attractor 

A.l Nearly-uniform distribution 

As a specific attractor consider the nearly-uniform distri- 
bution {pi{a)}^i: 



1 — a, i — 1] 
a/(iV- 1), i = 2..N. 



a G (0, 1) (70) 



The temperature ([8]) in the quasi-static limit ^ can 
be written as 



T = 5(7)eln 



_i (iV-l)(l-a) ^ 
a 



(71) 



with e := (iV - 1) ^Yln=2^i ~ ^i- From ei < it 
follows that e > 0, so the area of positive temperatures on 
the (Q;,7)-plane is restricted by a < 1 — N~^. 

The function 0, that defines the shape of isothermal 
paths 5(7) = T • 0{a), is now 



c'w(a) = e In : 



(72) 



and an example of the corresponding Carnot 's cycle is 
sketched in Fig. [11] 



Y 












Th ~ 








T 






i 
1 
1 


a 



ai 



Fig. 12. Carnot cycle in the (a, 7)-plane associated with Tsal- 
lis' distribution (I73|) and g{'y) = 7 at g = 1.4 and {e^} = 
{—3/2; —1/2; 1/2; 3/2}. The cycle is composed by segments of 
isentropes a = ai,2 and isotherms 7 = Tc,h • &Ts{a) (fTG]) ; 
(Ti > Tc). Note that the domain of definition for (|73[) is re- 
stricted by a < [maxi (q — . 



A. 2 Tsallis' distribution 

The attractor {pi{a)}fLi is now taken to be the Tsallis' 
distribution l28l: 



Pi{a) = Z^li^l - a{q - l)e, 



(73) 



with the specific parameter q, bounded by q € (0,2); at 
q = I the case reduces to the canonical one (fTS]). 

Introducing the mean value (•) and the covariance ((• ; •)) 
in a standard way: 



as well as an auxiliary variable 

a{q - l)ei 



Ui{a;q) 



1- a{q- l)€i' 



(74) 



(75) 



one can now express the isotherm shape function 0q-s (a; q), 
resulting from ([8|), as 



0rs{a;q) = a 



((u; In (1 + 7.))) 
{{u;u/{l + u)))' 



(76) 



Note that In(l-l-u) and u/(l -I- u) are both increasing 
functions of u, therefore the RHS is always positive here. 

An example of a Carnot cycle with isotherms of such 
a kind is sketched in Fig. [121 One can see that these 
isotherms are very similar to the canonical ones (cf . Fig. [T]) 
as long as a is small enough. This is easily shown taking 
aci < 1 for ([7S)1 . 



B Inequalities 

Here we prove some inequalities concerning the heat and 
work flows in case of the canonical attractor {pi{a)}fLi (H 
Considering the ffi'st derivative of Pi{a), one easily gets 



Xl^^'^ ~ Q)Pfc(a)Pi(a) 
da ^ 



(77) 
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and, consequently: with 

2—1 i,k—l 

„ x-^ , ,T / I \ One can easily check that both f ^ h(a)da = and 

CanZ^ltfc e ^u. (d/i/da) < hold. It follows from ^ that Y.^iV^{<^) 

is a decreasing function of a as well. The lemma proven 

So, one gets above ensures the integral in ([5^ to be positive and con- 

^e,(dp,/da) < 0, (78) sequently: 

that implies for a\ < the following inequalities: 



E Consider, at last, the term 'Y^tiXi from (f55l) . The in- 
(Pi(ai) - Pi(a2)) > 0; (79) ^^g^.^^ defined in ^ and reduces to: 



(Pi -Pi(a))ada = ( / + / ) (a - Q;)pi(Q!)da 



^e, / 'a dp, (a) < 0. (80) 

J OL\ 

Consider the terms V fpj^^') - p,(ai.2)) from m, 'T^''''' ^ = (^2 + ai)/2. It follows from the mean-value 

V ' / theorem that 

where the distributions at the cycle turning points _ 

are defined in m ■ Integration by parts leads to: ( J"'^ _ ^.p, (a)) da 

^ Jqi ^ ^ =^e^Pi[ei) / (a - a)da 

ra-2 

exp{(a - a,,,)/nr^} ~ exp{(a2;i - a)/A.r^.} + E ^MO^) j^^ {a - a)da 

exp{TZ\aAr^}-exp{±Z\«/A^r^,} ' = E (p.(0i) - ^,(^2)) ^ ^ d^. 

Observing Ti:2 > and ([75)l . one comes up with: 



where 



E/ (i;2) - . \\ ^ n /oi\ where ai < < a < 02 < a2. With one gets: 

(Pi -P»(ai;2)j ^0; (81) 

For the next inequality we need the following Aa ej^ = / (p^ — Pi(a)) a da > 0. (84) 



Lemma : For any two monotone decreasing functions J{x) 

rb 

an 



d h{x) the inequality h{x) f {x)dx > follows from q Asymptotes 

Jl^h{x)dx = 0. 

To get the asymptotic solution of the relaxation equa- 
Proof: Define H{x) /^^ h{x')dx' with the obvious prop- tion ([37]) in the slow driving limit {vTj^/Aa) < 1, we start 



erty: H{x) > H{a) = H{b) ~ for a < x < b. Applying an iterative procedure: For this purpose we rewrite ([3^ 
integration by parts, one gets: as 

p, = p,(a) - uTfl (dp,/da) ; (85) 

/ h{x)f{x)dx = [H{x)f{x)]l - / — if (a;)da; The zeroth approximation is given by pi « pi(a). The 

" b ° ^^^^ order is obtained by the corresponding substitution 

— rTt\A^ri right hand side of ((55|) : 



Pi~Pi- vtr {dpi /da) . (86) 



= - / ^i?(a;)dx > 
Ja dx 

observing {df /dx) < and i?(a;) > 0. 

This result, in turn, will be used, in the next iteration 
Consider now the term ^ (pf ^ — ) which occurs step: 

in RHS of (gH) and §7^/. With ^ it reduces to , o~ /j ^ . / ^2 (-,2- / - 2\ ^o^^ 

^. _ ^. _ ^^^^^^.^^^^ _^ ^^^^j (d Pi/da j , (87) 

e, (pf ^ - pf ^) = ^ |l - I (82) and so on. 

X f"'^ h( ) ~( )] d Consider now the asymptotes for the distributions p^^'^^ 

^'^i ^^P''^'^) J at the cycle's turning points first, in the limit of very 



Jan Birjukov et al.: Quantum thermodynamic processes: A control theory for machine cycles 



13 



slow driving, here [nTj^/Aa) ^ 1. Both p^P and pf''^ can 
be regarded as particular solutions of Eq. (ISS)) obtained at 

the specially chosen parameters and initial values. So, 
is achieved at a = ai through the evolution with v = —Xk 
and 
gets: 



and ao = ^2 ^ started at pf^ = pf \ With ^ one 



pf ^ = Pi(ai) + Akt^ 



(^) +o[ 



'R\ ' 



(2) 

Treating in the same way , one comes up with 



(88) 



(89) 



In the limit of very fast driving, {Aa/ kt^) <^ 1, one 
has to expand the exponents contained on the RHS of (|48l) 
into power series. Collecting terms up to first order in 

yields the results for p*^^'^^ presented in Eq. ((52|) . Their 
difference, (p^^-* — p*-^^), however, turns out to be of the 
second order and needs the adequate accuracy in the treat- 
ment. 



D Proof for Ql^ < 

In the stationary cyclic regime the system state evolution 
on the low-temperature stage 1 ^ 2 is given by Eq. p8|) 



taken with u = k, ao = ai and pf ^ = pl'^'' . Let us denote 



this particular solution of the relaxation equation (|37|) 
as p*; its boundary values prove to be p^(q;i;2) = p] 
from (|48)) . Consider the derivative of the sum J2 ^iP* ^ 
function (p{a): 



(p{a) 



da 



^e,p* -(kt^) i^e,(p* -p,(a)) 



applying 

It follows from ([8T|) that the boundary values of f{a) 
turn out to be (p(ai;2) ^ 0, so there must be at least one 
a* £ (ai;a2) such that >f{a*) = and, consequently, 

^e,p-(a*) ^ ^e,pi(a*). 

Choosing a* as a new initial point for (j38p . one can employ 
this equality in order to recast ^ e^p* and then </?(«) into 
the new form: 

^<^iP* = ^e^Piia*)e~^^ 

+ {i^TRy^ / w2<^iPi((^')) e'^^da' 
= ^<^iPi{oi)- I e"^^d f^e.^p,') ; 

J a* 

(^(a) = (kt^)"^ y e~^^d e,p^ . 

With ([71]) it is clear now that 

V3(a)^0 for a^a*, (90) 



so a* introduced above is unique. 

Consider now the integral originated from 

I "df^e.pn = / aip{a)da, 



which yields the total heat exchange on the low- 
temperature isotherm (7(7) — aTi of the Carnot cycle. 
It is convenient to split the interval of integration by a* 
so that the mean-value theorem becomes employable: 

+ I ) a {p{a) da — 9i / Lp{a)da 

+ 62 (p{a) da, 

J a* 

where ai < 61 < a* < 62 < a2. According to ([90| . the 
first addend here is positive while the second one is nega- 
tive. One estimates further: 



(^1 



^"2 

72- / ) ip{a) da 



< 92 / (^(a)da = 02 / d[Ve,p* 



as it follows from (|83p . 

This completes the proof for 



< 0. 



(91) 
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